Particle Ratios as a Probe of the QCD Critical Temperature 
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We show how the measured particle ratios can be used to provide non-trivial information about 
the critical temperature of the QCD phase transition. This is obtained by including the effects 
of highly massive Hagedorn resonances on statistical models, which are used to describe hadronic 
yields. The inclusion of Hagedorn states creates a dependence of the thermal fits on the Hagedorn 
temperature, Th, which is assumed to be equal to Tc, and leads to an overall improvement of thermal 
fits. We find that for Au+Au collisions at RHIC at ^/snn = 200 GeV the best square fit measure, 

, occurs at Tc ~ 176 MeV and produces a chemical freeze-out temperature of 172.6 MeV and a 
baryon chemical potential of 39.7 MeV. 



I. INTRODUCTION 

Lattice QCD is the main non-perturbative theoreti- 
cal tool used to probe bulk thermodynamics quantities 
of QCD such as its pressure, entropy density, and the 
speed of sound. The QCD phase transition at vanish- 
ing baryonic chemical potential is a (rapid) crossover 
where the thermodynamic quantities vary significantly 
near a critical temperature, whose value lies between 
170 - 200 MeV. In fact, according to the Bielefeld- 
BNL/RIKEN-Columbia collaboration (RBC-Bielefeld) 
the critical temperature is around Tc = 196 MeV [l| 
(although recently it has been concluded that the range 
could be Tc = 180 - 200 MeV 0) whereas the Bu- 
dapest/Marseille/ Wuppertal (BMW) collaboration has 
found a lower value Tc = 176 MeV Because the 
value of the critical temperature is vital to many phe- 
nomenological models of QCD, it is clearly important to 
find experimentally driven signals able to distinguish be- 
tween these two critical temperature regions. We shall 
show in this Letter that thermal fits for the measured 
particle ratios in Au-|-Au collisions at ^/snn = 200 GeV 
at RHIC can be used to determine the critical temper- 
ature of the QCD phase transition at nonzero baryonic 
chemical potential as long as effects from highly massive 
Hagedorn resonances are included. 

Thermal fits computed within statistical models are 
normally used to reproduce hadron yield ratios in heavy 
ion collisions d, d, @, 0, H, @- Thermal models com- 
puted at AGS, SIS, SPS, and RHIC energies can be 
used to construct a chemical freeze-out line in the QCD 
phase diagram [lo|. For Au-|-Au collisions at RHIC at 
y/SNN = 200 GeV, specifically, estimates for the chemi- 
cal freeze-out temperature and baryon chemical potential 
range from T^h = 155- 169 MeV and ^6 = 20 - 30 MeV 

Bii. 

Originally, it was thought that the chemical freeze- 
out temperature and the critical temperature coincided. 
However, recent lattice results indicate a higher critical 
temperature, which leads to a difference of AT = 7 — 45 
MeV between Tc and Tch- At SPS this difference was 



explained by allowing hadrons, especially exotic anti- 
baryonic states, to be "born" out of chemical equilib- 
rium at Tc and reach chemical equilibrium through multi- 
mesonic collisions [H, [l^ because chemical equilibration 
times of binary collisions are too long [l3|. At RHIC, 
multi-mesonic collisions are no longer adequate to ex- 
plain chemical equilibration times (isl . [l6j and this has 
led some to believe that hadrons are "born" in chemi- 
cal equilibrium [13, HI] ■ A way out of this scenario in- 
volving an overpopulation of pions and kaons has been 
suggested in [T^. Another solution that has provided 
very promising results is the inclusion of Hagedorn states, 
which are heavy resonances with an exponentially grow- 
ing mass spectrum that open up the phase space 
and help drive hadrons quickly into chemical equilib- 
rium [lil, m, m, |23|._ When a reaction of the form 
riTT ^ HS ^ rnr + XX is used where XX = pp, KK, 
or AA, hadrons are able to reach chemical equilibrium 
at about Tch ~ 160 MeV [2^ . [23] using various lattice 
critical temperatures [1^. Moreover, it was shown that 
the K/tt and {B + B)/'k ratios, where B — p + n, match 
RHIC data well 

Not only have Hagedorn states provided a mechanism 
for explaining the temperature difference between Tc and 
Tch, but they have also been used to find a low r]/s in 
the hadron gas phase [25|. which nears the string the- 
ory bound r]/s — l/(47r) [26|. Calculations of the trace 
anomaly including Hagedorn states also fits recent lat- 
tice results well and correctly describe the minimum of 

near the phase transition found on the lattice [25j . 
Furthermore, estimates for the bulk viscosity including 
Hagedorn states in the hadron gas phase indicate that 
Q/ s increases near Tc, which agrees with the general anal- 
ysis done in [27| . 

Although Hagedorn's idea of an exponentially grow- 
ing mass spectrum originated in the late 1960's, recent 
experimental results maintain an exponential mass spec- 
trum albeit with a higher Tc [1^. Moreover, it has been 
recently shown in [29] that a Hagedorn spectrum does 
appear in QCD with a large number of colors. Moreover, 
thoughts on observing Hagedorn states in experiments 
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are given in [30| and their usage as a thermostat in [3ll |. 
A possible method for describing the cross sections of 
Hagedorn states was derived in [33 |. 

Since Hagedorn states have been shown to affect the 
chemical equilibration times, thermodynamic properties, 
and transport coefficients of hadron resonance gases close 
to Tc it is natural to expect that they may also be relevant 
in the thermal description of particle ratios. Moreover, 
because Hagedorn states are dependent on the limiting 
Hagedorn temperature Th = Tc, a relationship between 
the chemical freeze-out temperature and the critical tem- 
perature can be found by including Hagedorn states in 
thermal fits. This uniquely gives us the ability to dis- 
tinguish between different critical temperature regions 
depending on the quality of the fit obtained using the 
statistical model. 



II. MODEL 

In this paper we use a grand-canonical model to de- 
scribe the particle densities from which we can calculate 
the corresponding ratios as described in detail in Q . We 
do not include any strangeness suppression factor or, in 
other words, we assume 7^ = f . In order to calculate the 
baryonic chemical potential fib and the strange chemical 
potential /is we use the following conservation relation 



corrections '33 







(1) 



which means that the total strangeness per baryon num- 
ber is held at zero. There rii is the density of the i*'' 
particle that has a corresponding baryon number Bi and 
strangeness Si. 

Hagedorn states are included in our hadron resonance 
gas model via the exponentially increasing density of 
states 
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which follows from Hagedorn's original idea that you have 
an exponentially growing mass spectrum that has a lim- 
iting temperature, Th- Close to Th the Hagedorn states 
become increasingly relevant and heavier resonances are 
"formed" the closer you get to Th- We use the particles 
from the particle date group up until Afg — 2 GeV and 
then we use Hagedorn states above 2 GeV. Additionally, 
A is the Hagedorn state "degeneracy", M is the maxi- 
mum mass, and m,. — 500 MeV, which is taken from 

In this paper we use two different scenarios regarding 
Th- First we assume that Th — Tc, and then we consider 
the two different different lattice results for Tc- Tc = 196 
MeV m and Tc = 176 MeV Q. Futhermore, we take 
into account effects from repulsive interactions between 
the hadrons [H, [s^l via the following excluded- volume 
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Note that the system's temperature T, baryonic chemical 
potential fib as well as the thermodynamic functions (af- 
ter volume corrections) are defined in terms of the quan- 
tities computed in the point particle (subscript pt) ap- 
proximation (i.e., no volume corrections). Note that B is 
equivalent to an effective MIT bag constant and is taken 
as a parameter in our model. 

In order to find the maximum M Hagedorn state 
masses and the degeneracy A, we fit our model to the 
thermodynamic properties of the lattice at zero chemi- 
cal potential /if, = 0. In the RBC-Bielefeld collabora- 
tion the thermodynamical properties are derived from 
e — 3p, which is what we fit in order to obtain the pa- 
rameters for the Hagedorn states. In this case we set 
Th = 196 MeV and A = 0.5 GeV^^^ and obtain M = 20 
GeV and B = (340 MeV)**. For the BMW collabora- 
tion the energy density is fitted and we fix Th — 176 
MeV and A = 0.5 GeV^^^ and obtain M = 15 GeV and 
B = (250 MeV) . Additionally, a discretized version of 
the resonance spectrum is considered, which is separated 
into mass bins of 100 MeV. Only mesonic, non-strange 
Hagedorn states are considered in our calculations. 

In our model we do not just consider the direct number 
of hadrons but also the indirect number that comes from 
other resonances. For example, for pions we consider also 
the contribution from resonances such as p's, w's etc. The 
number of indirect hadrons can be calculated from the 
branching ratios for each individual species in the particle 
data book Moreover, there is also a contribution 

from the Hagedorn states to the total number of pions, 
kaons, and so on as described in [13, H^. Thus the total 
number of "effective" pions can be described by 
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whereas the total number of "effective" protons, kaons, 
or lambdas (generalized as X) can be described by 



Nx = Nx 
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where {X) is the average number oi X = p's, K's, or A's. 
Here N is the total number of each species and (n^) is 
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the average number of pions that each Hagedorn state 
decays into. 

To determine {X) we use the multiplicities in Fig. 2 of 
Ref. [2lj from the microcanonical model in [s^ such that 



p = 0.058 mi -0.10 
K+ = 0.075 mi + 0.047 
A = 0.04 m,- - 0.07. 



(6) 



Clearly, they are all dependent on the mass of the i*^ 
Hagedorn state. Of course, in principle, the branching ra- 
tios of potential Hagedorn states are not known. Future 
measurements of high exotic hadronic resonances can be 
used to obtain these ratios in the future. Following the 
principle of (maximum) missing information, we assume 
here that the branching ratios can be obtained from a 
microcanonical calculation. Such a description is, for in- 
stance, also appropriate for describing the annihilation 
of p and anti-p. 

In order to get an idea of the quality of the thermal 
fits, we define a-s 
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where i?*''^'"™ is our ratio of hadron yields calculated 
within our thermal model whereas R^^^ is the exper- 
imentally measured value of the hadron yield with its 
corresponding error af. In this paper we look at only 
the experimental values at mid-rapidity and we used only 
the systematic error given by each respective experiment. 
We vary the temperature and according to the con- 
servation laws in Eq. ([T|), in order to get the smallest 
yf. We use the experimental data from both STAR 
and PHENIX W\ for Au-f Au collisions at RHIC 
at ^sjvAr=200 GeV. Specifically, we observe the ratios: 
tt^/tt+, p/p, K^/K+, K+/tt+, p/tt+, and (A + A)/tt+. 
All of which are calculated by STAR 37]. However, 
only TT"/?!""*", p/p, K~/K~^, K^/n^, p/n^ are given 
by PHENIX. Because there is such a difference between 
p/tt^ from PHENIX and STAR we choose only the value 
from STAR so that we can compare are results to [8] 
where they also exclude p/7r+ from PHENIX. It should 
be noted that Ref. @ includes more ratios than we do 
such as multi-strange particles and resonances, which are 
not included in this paper. This is because the purpose 
of this paper is not to confirm their results, which have 
already been confirmed in Q , but rather to compare ther- 
mal fits that include the contribution of Hagedorn states 
and those that exclude them. 



III. RESULTS 

The following results are given for the minimal foi' 
a given fit and Tch- Initially, we found the thermal fit for 
a hadron gas excluding Hagedorn states, which is shown 
in Fig. [TJ There Tch = 160.4 MeV, and ^f, = 22.9 MeV, 



which gave = 21.2. Our resulting temperature and 
baryonic chemical potential are almost identical to that 
in dl where T^h = 160.5 and fib = 20 McV. 

RHIC 200 GeV Au+Au 
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□ PHENIX 
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FIG. 1: Thermal fits for Au-I-Au collisions at RHIC at 
^/sNN = 200 GeV when no Hagedorn states are present. 

The inclusion of Hagedorn states is our primary inter- 
est. Starting with the fit for the RBC-Bielefeld collab- 
oration, we obtain Tch = 165.9 MeV, fit = 25.3 MeV, 
and — 20.9, which is shown in Fig. [2l The x^ is actu- 
ally slightly smaller than in Fig. [T] The contributions of 
the Hagedorn states to the total number of the various 
species at this temperature and chemical potential are 
shown in Tab. ID 

RHIC 200 GeV Au+Au 

with HS 



10 



10" 



10" 



10" 




+ T= 1 65 .9 MeV |i=25 .3 MeV Z"=20.9 (T_,= 1 96 MeV) 
T=172.6MeV |i=39.7 MeV x"=17.8 (T^-176MeV) 



7i"/7i"' p/p K /K"" K"'/7l'' p/7l"'(A+A)/7l" 



FIG. 2: Thermal fits including Hagedorn states for Au+Au 
collisions at RHIC at ^ilTiv = 200 GeV. 

When we consider the lattice results from BMW, which 
are at the lower end of the critical temperature spectrum 
where Tc = 176 MeV, we find = 172.6 MeV, fn, = 
39.7 MeV, and — 17.8. The lower critical temperature 
seems to have a significant impact on the thermal fit. The 
lower x^ is due to the larger contribution of Hagedorn 
states at at Tch = 172.6 MeV, which is much closer to 
Tc- The contributions of the Hagedorn states to the total 
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Te (MeV) 


A (GeV^/') 


M (GeV) 


tt's 


K's 


p's 


A's 


176 


0.5 


15 


48.5% 


41.6% 


29.1% 


41.0% 


196 


0.5 


20 


11.2% 


10.5% 


4.7% 


6.2% 


176 


0.5 


30 


48.7% 


41.6% 


29.1% 


41.2% 


196 


0.5 


40 


11.2% 


10.5% 


4.7% 


6.2% 


176 


1.0 


15 


62.5% 


56.0% 


40.6% 


53.4% 


196 


1.0 


20 


44.0% 


38.9% 


21.9% 


30.3% 



TABLE I: Contribution of the Hagedorn states to the total 
number of hadron species. 



number of the various species at this temperature and 
chemical potential are about 30 — 50% as shown in Tab. 
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RHIC 200 GeV Au+Au 

with HS 
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FIG. 4: Thermal fits including Hagedorn states for Au+Au 
colhsions at RHIC at ^/snn ~ 200 GeV when the degeneracy 
of the Hagedorn states is doubled. 



n (MeV) 


A (GeV^/^) 


M (GeV) 


(MeV) 


MS (MeV) 




176 


0.5 


15 


172.6 


39.7 


17.8 


196 


0.5 


20 


165.9 


25.3 


20.9 


176 


0.5 


30 


172.6 


39.8 


17.8 


196 


0.5 


40 


165.9 


25.0 


20.9 


176 


1.0 


15 


172.3 


51.7 


20.4 


196 


1.0 


20 


182.2 


33.5 


18.4 



TABLE 11: Comparison of the chemical freeze-out temper- 
ature, baryonic chemical potential, and for various fits 
including Hagedorn states. 



FIG. 3: Thermal fits including Hagedorn states for Au+Au 
collisions at RHIC at y^sjvjv ~ 200 GeV when the maximum 
mass of the Hagedorn states is doubled. 

The difference in the x^'s for BMW and RBC- 
Bielefeld collaboration is directly related to the contri- 
bution of Hagedorn states in the model. Because the 
RBC-Biclcfeld critical temperature region is significantly 
higher than its corresponding chemical freeze-out tem- 
perature the contribution of the Hagedorn states is min- 
imal at only 4-11% (see Tab. HI). To further prove this 
point we can vary the parameters that define the influ- 
ence of Hagedorn states in the model. If, for instance, 
we double the maximum mass we see in Fig. [3] that our 
thermal fits are not affected. This effect arises the true 
limiting temperature after volume corrections is larger 
than the critical temperature [1^ . The effects of chang- 
ing the maximum mass are only seen at temperatures 
larger than the critical temperature, which are not con- 
sidered in this study. 

While the maximum mass does not affect the quality 
of the fit, the parameter A, which is essentially the de- 
generacy of the Hagedorn states, does. The results of this 
are shown in Fig. [4] If we double A then we find that 
the minimum for Tc = 196 MeV has dropped down 
to = 18.4, which is only slightly higher than the best 
fit for Tc = 176 MeV in Fig. H This indicates that at 



Tc = 196 MeV more Hagedorn states would be needed 
in order to get a better fit. However, we also see that 
for Tc = 176 MeV and A=1.0 GeV^/^ that = 20.4. 
The reason for this is that there is an overpopulation of 
Hagedorn states. If we look at the contribution of Hage- 
dorn states to the individual particle species we see that 
the optimal contribution of Hagedorn states is around 
« 40 ± 10%, which is what we get for the fits = 176 
MeV, A=0.5 GeV3/2, M=15 GeV and Tc = 196 MeV, 
A=1.0 GeV3/2, M=20 GeV as seen in Tab. HI 

A summary of our results is seen in Tab. [Til We find 
that the inclusion of Hagedorn states should not only 
provide a better fit but they also affect the chemical 
freeze-out temperature and the baryonic chemical poten- 
tial. The more mesonic Hagedorn states are present the 
larger /if, becomes. Furthermore, our fits also have higher 
Teh's than seen in the fit without the effects of Hagedorn 
states. 



IV. CONCLUSIONS 

In this paper we assumed that the particle ratios mea- 
sured in Au+Au collisions at RHIC at ^/snn = 200 
GeV admit a purely statistical description at chemical 
freeze-out. Our results for thermal fits without Hagedorn 
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states concur well with other thermal fit models f§\ where 
the chemical freeze-out temperature {T^h = 160.4 MeV) 
is almost identical and the baryonic chemical potential 
{fit = 22.9 MeV) is only slightly larger. The thermal fit 
with the known particles in the particle data group pro- 
vides a decent fit with = 21.2. However, the inclusion 
of Hagedorn states provides an even better fit to the ex- 
perimental data. In fact, we find = 17.8, Tch = 172.6 
MeV, and fib = 39.7 MeV for the BMW collaboration 
while for the RBC-Bielefeld collaboration we obtained 
= 20.9, T^h = 165.9 MeV^ and j^j, = 20.9 MeV. This 
provides further evidence [l^, [H, [2J, [2^ that Hagedorn 
states should be included in a description of hadronic 
matter near Tc- Since the chemical freeze-out temper- 
ature was found to increase from 160 MeV to roughly 
165 MeV (RBC-Bielefeld) or 172 MeV (BMW) when in- 
cluding Hagedorn states, this exemplifies the degree of 
uncertainty in extracting chemical freeze-out thermody- 
namical parameters by means of such thermal analyzes. 

Furthermore, because Hagedorn states provide a 
bridge between the chemical freeze-out temperature and 
the critical temperature, we were able to make a quali- 
fied statement about which critical temperature region is 
more appropriate according to the quality of the thermal 
fits. We find that lower critical temperature regions are 
favored because more Hagedorn states are present close 
to the chemical freeze-out temperature and that a sub- 
stantial number of Hagedorn states (i.e. a contribution 
of about 40% to the total particle numbers) are needed 



in order to provide the best fit to the hadron yield ratios. 

A lower can be obtained for the higher critical tem- 
perature region when we double the degeneracy of the 
Hagedorn states, which would lead to a mismatch be- 
tween our thermodynamic quantities and those computed 
on the lattice (recall that the parameters that define the 
exponential spectrum in this case are obtained by fitting 
the results of the RBC-Bielefeld collaboration at = 0). 
As we can see from Tab. a change in the parameters 
even when they are doubled still gives a better fit than the 
thermal fits without Hagedorn states because a contribu- 
tion of Hagedorn states as small as 4-11% still contribute 
enough to lower x^- Therefore, this reconfirms the im- 
portance of including Hagedorn states in the hadron gas 
phase and, consequently, in the computation of thermal 
fits. Moreover, our results indicate that hadronization 
and chemical equilibration do not need to occur at the 
same temperature in order to explain RHIC data. 
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